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SECTION  II  -  INTERIM  REPORT 


Section  11  -  Interim  Rqport 

Method  of  probe  clouds  for  modeling  of  the  aerosol  dispersion  in  turbulence 
atmosphere  on  complicated  landscape 

This  interim  report  devoted  to  preparing  the  effective  methods  of  prediction  of  aerosol 
transport,  which  modifies  concentration  field.  This  method  takes  into  account  local 
convection  flows  and  whirlwinds  near  obstacles.  Movement  of  the  aerosol  plume  near  source 
depends  from  different  air  streams.  So  the  most  important  part  of  modelling  is  description  of 
aerosol  streams  and  their  route  from  the  source. 

Present  investigation  devoted  to  the  new  method  for  fast  and  accurate  direct  solving  of 
3D  aerodynamics  problem  with  complex  boundary  conditions.  The  synthetical  modelling  of 
aerosol  dispersion  in  atmosphere  is  based  upon  the  library  of  streams.  Library  contains  three 
dimensional  description  of  air  flows.  The  library  of  flows  must  be  big  enough  to  contain 
flows  around  various  kinds  of  obstacles  with  high  precision.  Each  flow  is  modelled  using 
final-differences  scheme  for  three  dimensional  Navier-Stokes  equations  (without 
temperature).  Final-differences  scheme  precision  could  be  improved  by  dividing  dynamics 
equations  by  value  of  average  density  and  introducing  the  kinematic  viscosity  of  air.  The  last 
replaces  theoretical  artificial  viscosity  so  it  is  not  needed.  Preparation  of  library  takes  a  lot  of 
time.  But  when  library  is  ready  computational  time  needed  for  synthetic  model  is  small 
enough  to  let  use  of  this  model  during  outdoor  experimental  studies. 

New  models  of  turbulence  diffusion  are  established. 

Orie  model  uses  the  streamline  which  passes  through  place  of  aerosol  generator  as  a 
center  line  of  Gaussian  model  of  aerosol  plume  dispersion.  This  modification  of  Gaussian 
model  is  used  for  evaluation  of  the  aerosol  dispersion  in  atmosphere  with  complicated 
orography.  Gaussian  models  provides  quasi  -  homogeneous  aerosol  distribution  inside 
aerosol  plume.  This  model  is  good  when  the  threshold  is  negligible.  But  it  can  not  be  used 
when  non  linear  effects  are  investigated. 

Model  of  another  type  is  based  on  the  method  of  probe  clouds.  Method  of  probe  clouds 
is  similar  to  method  of  markers  in  gasdynamics.  All  probe  clouds  come  from  nozzle  of 
aerosol  generator  and  move  along  the  streamline.  Every  cloud  have  it’s  own  parameters  for 
the  defined  set  of  physical  characteristics.  The  last  is  loaded  from  the  air  flows  library  for 
each  moment.  In  this  method  aerosol  plume  is  replaced  with  set  of  probe  clouds.  During 
simulation,  the  movement  of  each  cloud  is  divided  into  several  steps.  On  each  step  the 
corresponding  parameters  are  changed  to  model  elementary  physical  process.  For  example, 
the  stochastic  generation  of  clouds,  transformation  by  turbulence  diffusion  and  convectional 
transfer,  at  small  time  step  can  be  investigated.  The  global  field  of  any  physical  parameter  can 
be  found  from  description  of  probe  clouds  displacement  and  their  characteristics.  The  total 
number  of  clouds  can  be  chosen  to  reach  given  precision  or  to  compute  no  longer  the  jiven 
time  period.  This  method  is  effective  in  modelling  aerosol  transport  in  case  of  wheelwinds, 
convection,  etc.  This  investigation  enables  evaluation  of  complicated  problem  of  aerosol 
transport  during  field  experiment  in  quasi-real  time.  This  second  model  is  good  for 
computation  of  evolution  of  aerosol  with  given  precision.  It  is  very  important  for  aerosol 
optics,  investigation  of  other  non  linear  effects  in  real  atmosphere. 
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METHOD  OF  PROBE  PARTICLES  FOR  PC  MODELLING  OF  THE  AEROSOL 

DISPERSION  IN  TURBULENCE  ATMOSPHERE  FOR  COMPLICATED 

LANDSCAPE. 

•  N.N.  Belov,  P.N.  Belov 

Aerosol  Technology  LTD,  Russia,  fax  7-095-1474361,  Belov@Tehno.mmtel.msk.su 

KEYWORDS 

Atmosphere;  aerosol  dispersion;  streamlines;  3D  Navier-Stokes  equations 

PC  modelling  of  the  aerosol  plume  dispersion  in  atmosphere  is  very  important  for 
ecology.  It  is  very  important  to  find  real-time  information  about  transfer  of  dangerous 
aerosol  plume  in  atmosphere  after  incident  with  chemical  or  energetical  plant.  Results 
of  aerosol  field  experiments  closely  defined  from  information  about  aerosol  transfer  in 
atmosphere.  Great  number  of  human  activities  connected  with  aerosol  plumes,  with  its 
evolution  in  turbulence  atmosphere.  Most  important  and  most  difficult  for  evaluation  is 
zone  near  aerosol  source  (1-15  kilometres).  Inside  this  distance  aerosol  plume  closely 
depends  from  orography.  Inside  of  this  distance  aerosol  plume  has  small  dispersion.  It 
sizes  are  compatible  with  sizes  of  buildings  ,  mountains;...  Direct  calculation  of  this  gas 
dynamical  situation  is  extremely  complicated.  It  is  impossible  now  to  obtain  good 
estimation  of  this  task  by  finite  differences  method  PC  calculation.  So  most  important 
case  -  the  evaluation  of  the  aerosols  under  complicated  surface  and  near  aerosol 
generator  -  up  to  now  has  not  any  decision.  This  investigation  devoted  to  preparing 
effective  methods  of  investigation  of  transport  of  aerosol  plume  inside  of  near  -zone  (1- 
15  kilometres  from  source)  for  complicated  orography. 

Synthetical  model  of  aerosol  transport  modelling  for  the  complicated  orography. 

Most  difficult  problem  of  modelling  of  aerosol  transport  connected  with 
complicated  air  streams  near  buildings,  mountains  etc.  in  Boundary  layer.  The 
evaluation  of  aerosol  transport  may  be  calculated  by  several  methods  with  good 
precision,  Smirnov  (1996)  only  after  the  elevation  of  the  aerosol  clouds  by  convection 
and  large  whirlwinds  under  buildings  and  another  obstacles.  Moving  the  aerosol 
plume  near  source  depends  from  different  airstreams.  So  the  most  important  part  of  this 
modelling  is  way  of  aerosol  stream  between  and  near  buildings  ... 

Real  time  investigation  of  these  complicated  air  streams  is  impossible  before 
averawing  pf  the  aerosol  in  this  streamlines.  It  is  unknown  method  for  fast  and 
accurate  direct  solving  of  3D  aerodynamics  problem  with  complicated  boundary 
conditions. 

First  step  for  solving  of  this  problem  was  made  in  synthetical  modelling  of  aerosol 
dispersion  in  atmosphere,  Belov  (1998).  In  this  investigation  used  storage  of 
aerodynamics  stream  line  for  great  number  of  elementary  air  flows  for  different 
buildings,  group  of  buildings...  3D  computation  was  used  for  preparing  of  the  library 
of  the  flow  lines  for  great  number  of  obstacles.  This  library  prepared  BEFORE 
modelling  of  the  aerosol  plume  dispersion  provides  field  of  air  flows  for  great  number 
of  obstacles.  One  of  these  streamlines  which  moves  through  place  of  aerosol  generator 
provides  axes  line  for  Gaussian  model  of  aerosol  plume  dispersion. 

Well  known  Gaussian  models  for  evaluation  of  the  aerosol  dispersion  in 
atmosphere  are  fast  and  simple.  These  models  may  be  applied  for  great  numbers  of 
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types  of  atmosphere  conditions.  Wind,  convection,  temperature  distribution,  sun 
radiation...  -  great  number  of  parameters  may  be  used  by  this  model  with  good  results. 
But  all  calculations  may  be  fulfilled  for  plane  Earth  surface  only.  This  limitation 
reduces  of  significance  of  these  models.  Most  interesting  and  important  aerosol  plume 
problems  concern  with  people,  with  their  houses  and  industry.  These  places  suited 
mostly  inside  of  towns.  For  these  places  Gaussian  models  nobody  applied  before.  Our 
modification  of  the  Gaussian  model  provides  new  opportunities  for  this  class  of 
modelling  of  aerosol  dispersion  in  turbulence  atmosphere. 

Gaussian  models  provides  quasi  -  homogeneous  aerosol  distribution  inside  of 
aerosol  plume.  It  is  great  defect  of  this  model.  Most  important  properties  of  aerosol 
plume  (optical,  chemical...)  concerned  with  great  discontinuity  of  aerosols. 

Description  of  novel  aerosol  transport  modelling  for  the  complicated  orography. 

Most  important  results  of  our  investigation  concerned  with  spotted  technology  of 
description  of  the  evolution  of  aerosol  flow  in  turbulence  atmosphere.  Modelling  of 
the  evolution  of  the  aerosol  spots  is  very  important  for  aerosol  optics,  for  modelling  of 
aerosol  measurement  in  real  atmosphere.  For  this  purposes  it  was  used  method  of  probe 
clouds.  Each  spherical  small  cloud  generated  from  nozzle  of  aerosol  generation  moves 
by  streamlines  which  was  provides  from  streamlines  library.  Stochastic  generation  of 
these  clouds,  its  transformation  by  turbulence  diffusion,  convection  and  by  disturbing 
by  whirlwind  near  obstacles  was  used.  This  method  provides  good  results  with  standard 
Pentium  PC  and  CD  ROM  library. 

Conclusion 

This  investigations  provides  the  possibility  of  evaluation  of  more  complicated 
problem  of  aerosol  transport  in  turbulence  atmosphere.  Moreover  some  evaluation  of 
the  aerosol  dispersion  may  be  made  during  field  experiment  in  quasi-real  time. 
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Simulating  the  dispersion  of  aerosol  in  boundary  layer  of  atmospheric  is  a 
complex  problem  that  includes  analysis  of  various  air  flows  near  the  ground  and 
obstacles.  Route  of  aerosol  stream  strongly  depends  on  landscape  and  placement  of 
large  artificial  objects.  At  most  cases  Earth  shape  cannot  be  described  using  various 
coefficients  in  simple  equations  for  aerosol  transport.  Most  important  and  most 
difficult  for  evaluation  is  zone  around  .the  aerosol  source  in  range  of  2-15  kilometers. 
The  synthetic  model  of  aerosol  dispersion  [Belov(1998)]  in  atmosphere  simulates  aerosol 
transport  for  any  distance  if  all  air  flows  are  known.  Meteorological  measurements  are 
not  acceptable  for  this  model  because  of  too  few  details  of  three  dimensional  air  flows 
are  available.  Experimental  work  on  measuring  all  of  the  flows  velocity  and  structure  is 
hard.  To  produce  such  information  numerical  methods  are  used. 

The  present  method  is  based  upon  the  finite-differences  scheme  for  tree 
dimensional  Navier-Stokes  gasdynamics  equations  for  homogeneous  temperature  field. 
In  lower  atmosphere  absolute  values  of  variations  of  air  pressure  (and  density)  are  much 
less  than  the  average  value.  This  can  be  used  to  simplify  the  Navier-Stokes  equations: 
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Here  V=(v,u,w)  is  air  velocity,  p,^  -  average  value  of  pressure,  p^, 
average  value  of  pressure,  p  =  p„  +Pvr»  Piv  “  average  value  of  density, 
value  of  density,  p  =  p„+Pvr.il  -  dynamic  air  viscosity. 
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These  partial  derivative  equations  are  nonlinear.  To  solve  them,  to  find  air 
velocity  and  density,  final  differences  scheme  is  used.  In  order  to  reduce  scheme  round¬ 
off  errors,  in  equations  (l)-(3)  are  divided  by  density  variation  and  introducing  the 


air  kinematic  viscosity 


To  copmpute  values  on  each  new  time  step  the 


interation  process  is  initiated.  In  equations  (l)-(3)  nonlinear  terms  treated  as 
production  of  values  from  new  estimation  and  the  previous  one.  Estimations  are 
computed  by  solving  matrix  4x4  for  linear  algebraic  equations.  The  values  from 
previous  time  step  play  role  of  first  estimation. 

■  Shorted  form  for  the  scheme  is:  ■  = 

(5)  Vj  +  vv^  +uvy  +  vw,  =  -%•  (Pv,)s  +1 V-  +  v  (v5y  +  v^)+iv(u3,y  -h  w^-,) 

P  av 

(6)  Uj  +  VUj  -l-UUy  -t- wu^  =  (p„)y  +|v-  +  v(u5^  +u^)-i-iv-(v5J5^  -I- Wp) 


(7)  Wj  V  Wj  UWy  -I-  WW^  =  -  (Pvr)-z  +  4  V-  U^  +  V-  Wp)+  j  V  (v,^  +  U^) 

P  av 

(8)  (Pvr)t  +Pav  •  (Vj  +  Uy  +  W^)  +  V-  +  U*  (p„)y  +  W  (p„),  =  0 


Lower  index  x  mean  finite-differential  derivative.  Ex: 


V(i,  j,  k,  m)  —  V(i-  1,  j,  k,  m) 


This  method  was  successfully  implemented  and  tested  On  problem  of  horizontal 
wind  which  velocity  increases  with  height.  Bondary  conditions  can  be  of  three  kinds: 
1 -given  functions  (velocity  and  density),,  2-hard  surface,  3-free  air.  Boundary  conditions 
are  constructed  for  each  net  node  separate^ly.  Format  of  description  lets  add  new 
boundary  conditions.  Using  this  technique  non-usual  conditions  may  be  specified. 
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This  investigation  devoted  to  preparing  the  effective  methods  of  prediction  of 
aerosol  transport,  which  modifies  concentration  field.  This  method  takes  into  account 
local  convection  flows  and  whirlwinds  near  obstacles.  Movement  of  the  aerosol  plume 
near  source  depends  from  different  air  streams.  So  the  most  important  part  of 
modelling  is  description  of  aerosol  streams  and  their  route  from  the  source. 

Present  investigation  devoted  to  the  new  method  for  fast  and  accurate  direct 
solving  of  3D  aerodynamics  problem  with  complex  boundary  conditions.  The 
synthetical  modelling  of  aerosol  dispersion  in  atmosphere  is  based  upon  the  library  of 
streams.  Library  contains  three  dimensional  description  of  air  flows.  The  library  of 
flows  must  be  big  enough  to  Contain  flows  around  various  kinds  of  obstacles  with  high 
precision.  Each  flow  is  modelled  using  final-differences  scheme  for  three  dimensional 
Navier-Stokes  equations  (without  temperature).  Final-differences  scheme  precision 
could  be  improved  by  dividing  dynamics  equations  by  value  of  average  density  and 
introducing  the  kinematic  viscosity  of  air.  The  last  replaces  theoretical  artificial 
viscosity  so  it  is  not  needed.  Preparation  of  library  takes  a  lot  of  time.  But  when  library 
is  ready  computational  time  needed  for  synthetic  model  is  small  enough  to  let  use  of 
this  model  during  outdoor  experimental  studies. 

New  models  of  turbulence  diffusion  are  established. 

One  model  uses  the  streamline  which  passes  through  place  of  aerosol  generator  as 
a  center  line  of  Gaussian  model  of  aerosol  plume  dispersion.  This  modification  of 
Gaussian  model  is  used  for  evaluation  of  the  aerosol  dispersion  in  atmosphere  with 
complicated  orography.  Gaussian  models  provides  quasi  -  homogeneous  aerosol 
distribution  inside  aerosol  plume.  This  model  is  good  when  the  threshold  is  negligible. 
But  it  can  not  be  used  when  non  linear  effects  are  investigated. 

Model  of  another  type  is  based  on  the  method  of  probe  clouds.  Method  of  probe 
clouds  is  similar  to  method  of  markers  in  gasdynamics.  All  probe  clouds  come  from 
nozzle  of  aerosol  generator  and  move  along  the  streamline.  Every  cloud  have  it’s  own 
parameters  for  the  defined  set  of  physical  characteristics.  The  last  is  loaded  from  the  air 
flows  library  for  each  moment.  In  this  method  aerosol  plume  is  replaced  with  set  of 
probe  clouds.  During  simulation,  the  movement  of  each  cloud  is  divided  into  several 
steps.  On  each  step  the  corresponding  parameters  are  changed  to  model  elementary 
physical  process.  For  example,  the  stochastic  generation  of  clouds,  transformation  by 
turbulence  diffusion  and  convectional  transfer  at  small  time  step  can  be  investigated. 
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The  global  field  of  any  physical  parameter  can  be  found  from  description  of  probe 
clouds  displacement  and  their  characteristics.  The  total  number  of  clouds  can  be 
chosen  to  reach  given  precision  or  to  compute  no  longer  the  given  time  period.  This 
method  is  effective  in  modelling  aerosol  transport  in  case  of  wheelwinds,  convection, 
etc.  This  investigation  enables  evaluation  of  complicated  problem  of  aerosol  transport 
during  field  experiment  in  quasi-real  time.  This  second  model  is  good  for  computation 
of  evolution  of  aerosol  with  given  predsion;  It  is  veiy  important  for  aerosol  optics, 
investigation  of  other  non  linear  effects  in  real  atmosphere. 
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PC  modelling  of  the  aerosol  plume  dispersion  in  atmosphere  is  very  important. 
This  investigation  devoted  to  preparing  effective  methods  of  prediction  of  aerosol  . 
transport  in  area  around  the  aerosol  source  (1-15  kilometres)  for  complicated 
orography.  Known  models  of  aerosol  transport  can  not  be  realised  for  complicated  air 
streams  near  buildings,  excavations,  mountains  and  other  parts  of  landscape  has 
significant  influence  on  trajectory  of  air  streams.  On  small  distances  averaging  of  land 
surface  leads  to  serious  loose  of  precision.  Aerosol  transport  may  be  evaluated  by 
several  methods  with  good  precision  only  after  the  elevation  of  the  aerosol  clouds  and 
their  transported  over  the  obstacles.  Local  convection  flows  and  whirlwinds  near 
obstacles  modifies  concentration  field.  This  variation  is  too  small  if  the  source  of 
aerosol  is  far  from  the  investigated  area.  But  concentration  near  the  source  can  be  very 
high  at  one  place  and  there  can  be  no.  aerosol  at  the  another.  Movement  of  the  aerosol 
plume  near  source  depends  from  different  air  streams.  So  the  most  important  part  of 
modelling  is  description  of  aerosol  streams  and  thier  route  from  the  source. 

Present  investigation  devoted  to  the  new  method  for  fast  and  accurate  direct 
solving  of '  3D  aerodynamics  problem  with  complex  boundary  conditions.  The 
synthetical  modelling  described  in  Belov  et  all  (1998)  of  aerosol  dispersion  in 
atmosphere  is  based  upon  the  storage  of  descriptions  of  aerodynamics  streamlines  for  - 
great  number  of  elementary  air  flows  for  parts  of  landscape,  buildings  and  their  groups. 
The  library  of  streams  contains  three  dimensional  description  of  air  flows  for  great 
number  of  obstacles.  Each  flow  was  modelled  using  final-differences  scheme  for  three 
dimensional  Navier-Stokes  equations  (without  temperature)  which  is  described  in  Belov 
et  all  (1998).  Preparation  of  this  library  takes  a  lot  of  time.  But  when  library  is  ready 
computational  time  needed  for  synthetic  model  is  small  enough  to  let  use  of  this  model 
during  outdoor  experimental  studies. 

One  type  of  model  uses  the  streamline  which  passes  through  place  of  aerosol 
generator  as  a  center  line  of  Gaussian  model  of  aerosol  plume  dispersion.  This 
modification  of  Gaussian  model  is  used  for  evaluation  of  the  aerosol  dispersion  in 
atmosphere  with  complicated  orography.  Gaussian  models  provides  quasi  - 
homogeneous  aerosol  distribution  inside  aerosol  plume.  This  model  is  good  when  the 
threshold  is  neglibile.  But  it  can  not  be  used  when  non  linear  effects  are  investigasted. 


Behv  et  all  PC  modelling...  Contract  No:  681 7U97-M-571 2.  Page  2 


Model  of  another  type  is  based  on  method  of  probe  clouds  described  in  Belov  et 
all  (1998).  This  method  is  effective  in  modelling  aerosol  transport  in  case  of 
wheelwinds,  convection,  etc.  The  second  model  is  good  for  computation  of  evolution 
of  aerosol  cloud.  It  is  very  important  for  aerosol  optics,  investigation  of  other  non 
linear  effects  in  real  atmosphere.  Each  small  spherical  cloud  that  comes  from  nozzle  of 
aerosol  generator  moves  along  the  streamline.  The  last  is  loaded  from  the  air  flows 
library  for  each  moment.  Stochastic  generation  of  these  clouds,  their  transformation  by 
turbulence  diffusion,  convection  and  whirlwinds  can  be  easily  implemented.  This 
method  provides  impressive  results  on  Pentium  PC  and  library  on  CD  ROM.  This 
investigation  enables  evaluation  of  complicated  problem  of  aerosol  transport  during 
field  experiment  in  quasi-real  time. 
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